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We report that the ionic-hquid gating of bulk single crystals of a topological insulator can control 
the type of the surface carriers and even results in ambipolar transport. This was made possible by 
the use of a highly bulk-insulating BiSbTeSe2 system where the chemical potential is located close 
to both the surface Dirac point and the middle of the bulk band gap. Thanks to the use of ionic 
liquid, the control of the surface chemical potential by gating was possible on the whole surface of a 
bulk three-dimensional sample, opening new experimental opportunities for topological insulators. 
In addition, our data suggest the existence of a nearly reversible electrochemical reaction that causes 
bulk carrier doping into the crystal during the ionic- liquid gating process. 

PACS numbers: 73.25.-|~i, 73. 20. At, 72.20.My, 73.30.-|-y 



I. INTRODUCTION 

Three-dimensional (3D) topological insulators (TIs) 
are characterized by a novel topological order— ^— which 
dictates the appearance of spin-filtered massless Dirac 
fcrmions on the surface^"— To experimentally address the 
peculiar physics of 3D TIs, it is desirable to access the 
Dirac point of the surface state (SS)i^"— This is relatively 
easy with the surface-sensitive spectroscopies such as the 
angle-resolved photoemission^rJ^ and the scanning tun- 
neling microscopcr^^ but it is more challenging for the 
bulk-sensitive transport experiments because the chemi- 
cal potential is always pinned to the bulk bands (includ- 
ing the impurity band) in real materials. To tune the 
chemical potential to a desirable position for transport 
experiments, two approaches have been employed: one is 
the tuning of the chemical compositions upon synthesiz- 
ing crystalSfi^Ti^^ and the other is the gating to control 
the surface carriersi^"— Among the latter approach, the 
electric-double-layer gating (EDLG) method is a promis- 
ing new tcchnique^^ to allow application of a large electric 
fields 

In the EDLG configuration, either cations or anions in 
a liquid electrolyte are accumulated near the surface of a 
sample by application of an electric field, and they form 
an electric double layer which generates a very strong 
electric field locally on the surface. This leads to the in- 
duced surface carrier density of as high as ~10^^ cm~^i^ 
Such a large tunability of the surface carrier density 
is the merit of the EDLG method. In the context of 
TIs, Yuan et al. showed™ that the EDLG method al- 
lows an ambipolar doping control in ultrathin films of 
Bi2Te3; however, the surface state was obviously gapped 
in the ultrathin films used in Ref. [2^, and consequently 



the transport properties were likely to be dominated by 
the bulk state. In the present paper, we show that the 
EDLG method can achieve ambipolar transport even in 
bulk single crystals, if one uses the highly bulk-insulating 
Bis-^Sb^Tes-yScj; (BSTS) system-^i^i In our experi- 
ment, the control of the chemical potential was possi- 
ble on the whole surface of a bulk 3D sample, opening 
new experimental opportunities for TIs. In addition, our 
data suggest that the chemical potential is moving not 
only at the surface but also in the bulk, possibly due 
to a nearly reversible electrochemical reaction to cause 
apparent bulk doping during the EDLG process. 



II. SAMPLE PREPARATIONS AND 
MEASUREMENTS 

A series of BSTS single crystals were grown by a modi- 
fied Bridgman mcthodj ^°'^^ Six gold wires were attached 
to each sample by a spot welding technique, and mag- 
netotransport measurements were performed by a con- 
ventional ac six-probe method by sweeping the mag- 
netic field between ±9 T. An electrically insulating cup 
made of Stycast 1266 was used as a sample container, 
in which the sample with gold wires was submerged into 
an ionic liquid (IL) electrolyte, as schematically shown 
in Fig. 1(a). We used a specially purified ionic liq- 
uid [EtMelm] [BF4] as the electrolytcj^ As the gate elec- 
trode, an additional piece of gold wire was dipped into 
the IL without touching the sample. In this paper, we 
adopt the convention of defining the gate voltage by tak- 
ing the sample as the reference point, as was done in most 
of the previous works. In this convention, n-type carriers 
are supposed to be doped to the surface when a positive 
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FIG. 1: (color online) (a) Schematic picture of the experimental setup. Actual ionic liquid has no color, (b)-(e) Changes in 
parameters with time in the procedure of applying the gate voltage Vg with the target value of —3.5 V on a BSTS sample. 
First, Vg is gradually changed from to —3.5 V while keeping the temperature at 220 K, and we wait for at least 10 min for 
the electric double layer to develop; the system is then cooled down, and the gate current Ig vanishes completely after the ionic 
liquid solidifies, (f)-(h) Changes in T, Ig, and the integrated charge in the relaxation process from Vg ~ —3.5 V, in which Vg 
is set to V and the system is warmed to ~240 K. (i) Vg vs Ig curve at 220 K with the voltage sweep rate of ~— 0.5 V/min. 



gate voltage is applied. 

The gate voltage was applied in the following proce- 
dure: first, the temperature T of the sample was sta- 
bilized at 220 K; then, the gate voltage Vq was swept 
slowly, and the temperature was kept at 220 K for at least 
10 min after the voltage reached the set value; finally, 
the sample was cooled down slowly. When changing the 
gate voltage, we cycled the above procedure. Figures 
l(b)-l(e) show an example of the history of T, Vg, gate 
current Iq, and the accumulated charge, during the pro- 
cedure for the target Vg of —3.5 V. With the IL used in 
the present experiment, gate current was never detected 
at temperatures below 200 K, even when the gate voltage 
was changed. This is because when the IL solidifies, ions 
are not mobile at all. 

The accumulated ions can be released when the gate 
voltage is set to zero and the system is warmed up to 
220 K. Such a relaxation process from Vg = —3.5 V is 
shown in Figs. 1(f)- 1(h). One can see that the amount 
of released charge [~32 fj,C in Fig. 1(h)] is comparable 
to that of the accumulated charge [~34 /xC in Fig. 1(e)]. 
The small difference between the released and accumu- 
lated charge is probably an indication of an irreversible 
electrochemical reaction during the gating processes at a 
high gate voltage. In the above example, the total sur- 



face area of the sample was 13.4 mm^, so that the '^32 
/iC of charge accumulated on the surface corresponds to 
the accumulated ion density of ~1.6xl0^^ cm~^ and the 
capacitance of the unit area 80 /zF cm~^, which is com- 
parable to that previously reported for EDLGi^ Figure 
l(i) shows Vg vs Ig curve at 220 K, which indicates that 
there is a threshold |Vg| of '^2 V below which little cur- 
rent fiows, suggesting that the ions are not mobile below 
this threshold voltage; similar behavior was previously 
reported for a different ILjH The origin of this behavior 
is currently not clear, but it may well be a characteris- 
tic of the ionic liquid used here. Also, it was difficult 
to obtain reproducible results of the transport proper- 
ties for |Vg| between 2 and 3 V, probably because the 
formation of the electric double layer is unstable in this 
gate-voltage region. Therefore, we closely measured the 
transport properties only in the range of |Vg| > 3 V and 
Vg = V. 



III. EDLG EXPERIMENT ON Bii sSbo sTei rSei 3 
SINGLE CRYSTALS 



First, we show the results of the EDLG experiments 
on Bi1.5Sbo.5Te1.7Se1. 3, which is currently one of the 
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FIG. 2: (color online) (a),(b) Temperature dependences of 
the resistivity of a Bi1.5Sbo.5Te1.7Se1. 3 single crystal for Vc ~ 
—3.5, 0, and +3.5 V. (c) Magnetic-filed dependences of pyx 
at 1.8 K. In these three figures (a)-(c), the dotted, solid, and 
broken lines correspond to Vg of —3.5, 0, and +3.5 V, respec- 
tively, (d) Fitting of the two-band conduction model to the 
pyxiH) data at Vg = V, which yields pbuik = 3.2 SI cm, Tibuik 
= 9.1 X 10" cm-^ /ibuik = 29 cmVVs, p™ = 8.1 kfi, nss 
= 1.6 X 10^^ cm"^, and /iss = 4.7 x 10^ cm^/Vs. The gray 
thick line represents the experimental result and the solid line 
is the fitted result. 



most bulk-insulating TI materials 1^1^ Figure 2(a) shows 
the temperature dependences of the resistivity p with 
the gate voltage (Vg) of —3.5, 0, and +3.5 V. Even 
though this is a bulk single crystal, at low tempera- 
tures p presents a clear change upon application of Vg; 
namely, below ^50 K, p changes with changing Vg, and 
it tends to decrease (increase) for positive (negative) Vg- 
Figure 2(b) magnifies this change for T < 20 K. Fig- 
ure 2(c) shows the magnetic-field dependences of the 
Hall resistivity, pyx{H). The sign of the charge car- 
riers remains negative in this sample, although a sig- 
nificant non-linearity suggests the coexistence of sur- 
face and bulk conduction channels. Note that the sur- 
face conduction in TIs can involve both the topological 
surface states and topologically-trivial two-dimensional 
electron-gas (2DEG) states that appear as a result of 
band bending j^ In the present experiment, when the 
gate voltage is applied, the 2DEG states are most likely 
contributing to the conduction alongside of the topologi- 
cal surface states. Unfortunately, we have not been able 
to elucidate the contributions of the topological and non- 
topological surface states because no Shubnikov-de Haas 
(SdH) oscillation has been observed in our gated sam- 
ples. (To successfully separate the contributions of the 
two, detailed information obtained from SdH oscillations 
is necessaryi^SiM) We therefore make no claim of the com- 



position of the surface carriers in the present paper. 

In Fig. 2(c), one can see a clear tendency that both the 
absolute value of Py^ and the slope of Py^iH) decrease 
upon increasing Vg from —3.5 to +3.5 V; this means 
that the apparent electron concentration increases with 
increasing Vq- Since n-type carriers are expected to be 
doped when a positive Vg is applied, the present obser- 
vation can be understood as a natural consequence of 
the EDLG. The nonlinear //-dependence observed in py^ 
is useful for gaining insights into the respective roles of 
bulk and surface transport channels, because its analysis 
based on a simple two-band modeli^ gives a crude idea 
about the relevant transport parameters. For example, 
the Pyx{H) data at Vg = V give the following estimate 
based on the fitting shown in Fig. 2(d): For the bulk 
channel, the bulk resistivity pbuik — 3 f2 cm, the bulk 
carrier density nbuik — 9 x 10^^ cm""^, and the bulk mo- 
bility /ibuik — 30 cm^/Vs; for the surface channel, the 
sheet resistance pgg ~ 8 kfl, the surface carrier density 
ngs ~ 2 X 10^^ cm~^, and the surface mobility pss — 
5 X 10"^ cm^/Vs. In this fitting, the constraint imposed 
by the presence of sharp kinks at '-^ ±1 T helps reduce 
the ambiguity in the fitting parameters, and in fact, the 
three parameters, nbuik, Mbuik, and nss, are in reason- 
able agreement with our previous transport studies of 
BSTS involving SdH oscillations 1^2*21 The large value of 
/iss would imply that the SdH oscillations be observed, 
but we did not observe any SdH oscillations in this sam- 
ple; this is possibly because the surface chemical poten- 
tial (and hence nss) is not very uniform throughout the 
sample and the SdH oscillations with various different 
frequencies add up to smear visible oscillations. 

The above result of the two-band analysis suggests 
that the surface contribution in the total conductance 
was only '^1%, which is reasonable because the mea- 
sured sample had a considerable thickness of 332 pm. 
Nevertheless, if one looks at the EDLG effect in resis- 
tivity [Figs. 2(a) and 2(b)], one notices a very puzzling 
fact: for the negative Vg of —3.5 V, the number of sur- 
face electrons are expected to be reduced and indeed, the 
slope of Pyx gets larger; however, the resistivity increase 
is as much as 25%. Since the surface contribution in the 
total conductance is only ^1%, even when the surface 
conduction is completely suppressed by EDLG, one can 
expect an increase in p of '^1% at most, as long as the 
bulk channel is not affected by EDLG. Therefore, the 
observed large increase in p strongly suggests that the 
bulk channel must also be affected by EDLG. Indeed, the 
Pyx{H) data shown in Fig. 2(c) presents a clear change 
in the slope at high fields for different Vg, which implies 
that not only nss but also ribuik is changing. To corrobo- 
rate this inference, the two-band analyses of the pyx{H) 
data at finite Vg suggests that the bulk carriers decreases 
(increases) by 5% (10%) for Vq of -3.5 V (+3.5 V). 

It is very surprising that the EDLG affects the density 
of bulk carriers by a noticeable amount in a sample as 
thick as 332 /xm, but our transport data can hardly be 
understood if one does not accept this possibility. Given 
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TABLE I: Experimental order of applied Vq and resultant 
transport properties at low temperatures, 
field limit of the Hall coefficient. 



that the electric field generated by ~'l x 10^^ cm^^ of 
ions on the surface is shielded in less than 100 nm;^ the 
only possibility is that some bulk doping into the BSTS 
sample is taking place during the EDLG process. In this 
respect, the slow time scale of the change in the measured 
current during the EDLG process seems to support the 
idea that some electrochemical reaction is taking place. 
In passing, we note that we have not successfully mea- 
sured SdH oscillations in any of the gated samples. This 
is likely to be due to an inhomogeneous distribution of the 
local electric field (which is conceivable because our sam- 
ples have a lot of macroscopic terraces on the surface) or 
some chemical degradation of the surface caused by the 
ionic-liquid gating. Hence a more definitive analysis of 
the transport data beyond the simple two-band analysis 
is currently unavailable. Nevertheless, the bulk doping 
due to the EDLG seems to be an inevitable conclusion of 
our result. 



IV. AMBIPOLAR TRANSPORT IN BiSbTeSea 
SINGLE CRYSTALS 

It was recently foundi^ that the position of the chemi- 
cal potential with respect to the Dirac point of the SS is 
tunable within the bulk band gap in BSTS when one fol- 
lows a particular series of x and y that were identified in 
Ref. [2l|. In particular, the chemical potential was found 
to be close to both the Dirac point of the SS and the 
middle of the bulk band gap in BiSbTeSe2. Therefore, 
for our EDLG experiment, to maximize the possibility 
of achieving ambipolar transport, we mainly measured 
the BiSbTeSe2 system. (The gating data shown in Fig. 
1 were taken during the experiment on the BiSbTeSe2 
sample reported below.) 

Figure 3(a) shows the temperature dependence of the 
resistivity p of a BiSbTeSe2 single crystal before applying 
the gate voltage. One can see that p exceeds 10 Hem 
at low temperature, testifying to a high quality of this 
sample. As was the case with Bi1.5Sbo.5Te1.7Se1. 3, at 
low temperatures we observed clear change in p upon 
application of Vg [Figs. 3(b) and 3(c)] even though this 
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FIG. 3: (color online) (a) Temperature dependence of p in the 
BiSbTeSe2 single crystal used for the EDLG experiment at Vg 
= V in a semi-log plot; inset shows the Vg dependence of p 
at 1.8 K, where the arrows indicate the order of experiments. 
(b),(c) Low-temperature part of p{T) for various Vg- 



is a bulk single crystal. However, it also turned out that 
the transport properties do not completely recover after 
cycling Vg ■ To illustrate the situation, we show in Table 
I the measured transport properties at 1.8 K for various 
values of Vg in the order of the measurements. As one 
can see in this table, the transport properties for Vg = 
were measured twice, before and after applying Vg = 
—3.5 V; also, there arc two different data for Vg — —3.3 
V, which were taken before and after Vg — -1-3.0 V was 
applied. In both cases, the resistivity decreased after a 
high voltage was applied. Such a decrease in resistivity 
after application of a high voltage was also observed in 
other samples, so it appears to be an unavoidable effect in 
BSTS crystals; this is probably due to some irreversible 
electrochemical reaction taking place in the bulk, which 
gradually spoils the bulk-insulating property. However, 
as one can see in the inset of Fig. 3(a), a large part of 
the change in resistivity in response to Vg is reversible. 
This reversible part of the change is consistent with the 
picture that the number of n-type carriers increases with 
positive Vg and decreases with negative Vg due to the 
EDLG effect. 

An interesting feature in our resistivity data is that 
the resistivity value presents a maximum around —3.2 
V [see inset of Fig. 3(a) and also Fig. 3(b)]; namely, 
the resistivity is smaller at Vg of —3.5 V compared to 
that at —3.3 V, despite the overall trend that negative 
voltage increases p. In fact, when the system is in the 
n-type regime, a more negative Vg value would lead to a 
smaller number of n-type carriers, and one would expect 
the resistivity to increase; the opposite behavior observed 
for Vg < —3.3 V suggests that the system is changing 
from n-type to p-type. To confirm this possibility, one 
must look at the Hall data. 

Figure 4(a) shows how the behavior oi pyx{H) changes 
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FIG. 4: (color online) (a) Magnetic-field dependences of pyx 
at 1.8 K for various Va- (b) Vq dependence of Rh at 1.8 K. 
(c) Temperature dependences of Rh for various Va- Arrows 
indicate the experimental order. 



with Vg- Initially, at Vg = V (first) and at -3.125 
V; the charge carriers are clearly n-typc, but the carriers 
become p-type at large negative Vg values of —3.5 and 
—3.3 V. This is a signature of an ambipolar transport in 
a bulk crystal, and confirms the idea that the sign change 
of charge carriers takes place between Vg = —3.125 and 
—3.3 V. As noted above, the apparent maximum in the 
resistivity near Vg — —3.2 V is consistent with this in- 
terpretation. The n-type doping is recovered when Vg 
was set to V again, and the slope of pyx {H) was found 
to decrease with increasing Vg up to -1-3.0 V, suggesting 
an increase in the n-type carriers, as expected. 

It should be noted that we observed the sign change 
in this sample again after setting Vg to -1-3.0 V and then 
bringing it back to —3.3 V, as shown by a broken line 
in Fig. 4(a). Thus, the sign change of the carriers is 
obviously reproducible. By defining the Hall coefficient 
Rh as the slope of pyx{H) at low field, we summarize the 
gate- voltage dependences of Rh in Fig. 4(b). One can 
see that the sign change in Rh is reproducibly observed, 
although the exact value of Rh at a given Vg shows a 
history dependence. 

The temperature dependences of Rh for various Vg 
measured in the successive five experiments [Fig. 4(c)] in- 
dicate that the ambipolar transport is observed only be- 
low ^30 K where thermal activations of bulk carriers are 
negligible. Since the experiment on Bi1.5Sbo.5Te1. 7861. 3 
discussed in the previous section indicated that the car- 
rier densities in both the bulk and surface transport chan- 
nels are changing with EDLG, it is important to eluci- 
date whether the sign change of the carriers observed in 
BiSbTeSe2 is occurring in the bulk or on the surface, or 
both. To infer the origin of the sign change, we have an- 
alyzed the Pyx (H) data for various values of Vg with the 
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FIG. 5: (color online) Results of the two-band model fitting 
to the pyx{H) data of BiSbTeSe2 at various Vc- Thick lines 
are the data at 1.8 K and solid lines are the fitting results. 
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TABLE II: Parameters of the two-band model fitting of the 
pyx{H) data at various Vg; the fitted curves are shown in Fig. 
5. 



simple two-band model. 

The results of the two-band fitting are summarized in 
Fig. 5 and Table II. For example, the fitting to the data 
at Vg = V shown in Fig. 5 gives the following crude 
estimate for the transport parameters: For the bulk chan- 
nel, pbuik — 13 ricm, ribuik ^ 1 X 10^'' cm"^, and /Xbuik — 
40 cm^/Vs; for the surface channel, pgg — pss/d — 30 
kn, ngs ~ 2 X 10" cm^^ and ^gg ~ 1 x 10^ cm^/Vs. 
Here, d (= 181 pm) is the thickness of the sample and 
the sign of charge carriers is negative for both the bulk 
and surface channels. Given that the chemical potential 
is very close to the Dirac point in BiSbTeSe2, it seems 
that the estimate of the surface carrier density indicated 
in this analysis is reasonable, in spite of the weakness 
of the non-linearity in pyx{H). With the above parame- 
ters, the contribution of the surface channel to the total 
conductance is calculated to be ~2%. 

Looking at the results of the two-band analyses of the 
data at finite Vg, one can see that the Hall data strongly 
suggests that both the bulk and surface carriers change 
sign simultaneously at Vg < —3.3 V. Indeed, we found 
that it is impossible to fit the pyx data for Vg < —3.3 V 
by assuming that only one of the two channels changes 
sign. Therefore, it appears that in the present EDLG 
experiment the chemical potential is swung from the n- 
type regime to the p-type regime not only on the surface 
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FIG. 6: (color online) Magnetoresistance data at various 
stages of the EDLG experiment of BiSbTeSe2 at 1.7 K. 



but also in the bulk. Most likely, what is happening in 
the bulk at negative Vg is compensation due to electro- 
chemical p- type doping, which eventually overwhelms the 
preexisting n-type bulk carriers. 

Finally, we show in Fig. 6 the magnetoresistance 
(MR) data for various Vq values. Obviously, the data 
present the weak antilocalization (WAL) effect at low 



fields 



,34.39.40 



since the WAL effect is a signature of two- 



dimensional transport and is not usually observed in bulk 
TI crystals dominated three-dimensional transport, the 
MR data give additional evidence for a sizable contribu- 
tion of the surface transport to the total conductance. 
At the same time, this WAL effect makes it difficult to 
examine the consistency of the parameters obtained from 
the two- band analysis of pyx{H) in the MR data. Also, 
the behavior of MR qualitatively changed after the ap- 
plication of Vg = +3 V, the origin of which is not clear 
at the moment. Since the MR is too complicated and 
not very reproducible, we did not try to make a detailed 
analysis. 



V. DISCUSSION 

It is important to mention that there are missing 
charges in our EDLG experiment; namely, the total 
amount of charge induced by gating in the sample is much 
smaller than the total amount of ions accumulated on the 
surface. For example, in the case of BiSbTcSe2, the total 
amount of accumulated charges measured by the current 
is -^34 /xC, which corresponds to the ion density on the 
surface of -^1.6 x 10^'' cm^^. On the other hand, the 
change in the surface carrier density in BiSbTeSe2 was 
~4 X 10^^ cm~^ and its bulk carrier density changed 
by less than 1 x 10^^ cm^'^, which amounts to the to- 
tal charge of less than 20 ^G. Therefore, obviously the 
gating is not very efficiently performed. In this regard, 
the doping control of the surface carriers in the present 
experiment is similar to another EDLG experiment on 



Bi2Te3 thin film^^ where ^7 x 10^^ cm~^ of surface car- 
riers were doped with Vg — —3 V. For other materials, 
the amount of surface carrier doping by EDLG is of the 
order of lO^-lO^^ ^^-2 29^^i^^&^ ^^^^ ^j^^g ^j^g electric- 
field effect on Bi-based topological insulators appears to 
be exceptionally ineffectivei^ 

The mechanism to cause these missing charges is not 
clear at the moment, but we speculate that some elec- 
trochemical redox reaction involving adsorbed molecules 
on the surface causes a layer of immobile charges that 
shields some fraction of the electric field created by the 
ions, leading to a weakening of the electric field for in- 
ducing mobile carriers in the sample. 

Also, the bulk doping that accompanies the EDLG in 
our BSTS samples is surprising. Remarkably, the data for 
BiSbTeSe2 suggests that the bulk doping process is nearly 
reversible and it takes place in the time scale of the order 
of 10 min. The chemical mechanism of this bulk doping 
is not clear at the moment, but the possible cause might 
be intercalation of ions into the van-der-Waals gap in the 
BSTS crystal. Obviously, there is a lot to understand 
about the electrochemistry accompanying the EDLG on 
Bi-based tetradymite TI materials. 



VI. CONCLUSION 

In conclusion, the electric-double-layer gating (EDLG) 
using ionic liquid was applied to bulk single crystals of 
BSTS to control the chemical potential, and ambipo- 
lar transport was observed in a sample of BiSbTeSe2 
as thick as 181 /im. The gating was successfully ap- 
plied to tune the chemical potential on the whole sur- 
face of a three-dimensional sample, and surprisingly, it 
appears that the EDLG on BSTS crystals is accompa- 
nied by a nearly reversible electrochemical reaction that 
caused bulk carrier doping. It turned out that the EDLG 
is exceptionally inefficient for the BSTS system, with the 
maximum change in the surface carrier density of ^4 
X 10^^ cm~^ despite the ion density on the surface of 
~1.6 X 10^^ cm~^. The key to the successful ambipolar 
carrier control in the present experiment was the use of 
BiSbTeSe2 crystals in which the chemical potential is lo- 
cated close to the middle of the bulk band gai>i^ and the 
residual bulk carrier density was only ^1 x 10^^ cm"'^. 
In combination with a technique to open a gap on the 
surfacej^i^ the present experiment paves the way for 
topological magnetoelectric-effect experiments^ that re- 
quire the chemical-potential control on the whole surface 
of a bulk topological insulator, although the mechanism 
of the bulk doping associated with EDLG needs to be 
understood before this technique is comfortably applied. 
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